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NOMENCLATURE 


distance  in  semi -chord  measured  from  mid-chord  to 
elastic  axis 

matrix  of  aerodynamic  coefficients 
semi -chord  and  full -chord,  respectively 
critical  pressure  coefficient 
lift  coefficient 

moment  coefficient  about  pitching  axis 

moment  coefficient  about  aileron  hinge  axis 

distance  in  semi-chcrd  measured  from  mid-chord  to  aileron 
hinge  axis 

differential  c  relative  to  the  mean  value  of  c.  in 
forced  motion  * 

differential  c  relative  to  the  mean  value  of  cm  in 
forced  motion  m 

differential  c  relative  to  the  mean  value  of  c  in 
forced  motion  n 

structural  damping  coefficient 

plunging  degree  of  freedom 

polar  moment  of  inertia  about  elastic  axis 

polar  moment  of  inertia  about  aileron  hinge  axis 

wc/U,  reduced  frequency  with  respect  to  full  chord 

matrix  of  stiffness  coefficients 

aerodynamic  load  vector 

displacement  vector 
2  1/0 

(I^/mb  )  ,  radius  of  gyration  about  elastic  axis 

2  1/2 

(Ig/mb  )  ,  radius  of  gyration  of  aileron  about 

aileron  hinge  axis 

airfoil  static  moment  about  elastic  axis 
aileron  static  moment  about  aileron  hinge  axis 
time  in  seconds 
cat,  nondimensional  time 
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U  -  free  stream  velocity 

U*  -  U/bco^,  nondimensional  flight  speed 

x  S  /mb,  distance  in  semi-chord  measured  from  elastic  axis 

a  a 

to  mass  center 

Xg  -  S^/mb,  distance  in  semi-chord  measured  from  aileron  hinge 
axis  to  aileron  mass  center 

a  -  pitching  degree  of  freedom 

3  -  aileron  pitching  degree  of  freedom 

6  -  h/c,  nondimensional  plunging  degree  of  freedom 

<k  -  phase  difference  between  plunging  and  pitching  oscillation, 

’  plunging  leading  the  pitching 

4>h  -  phase  difference  between  plunging  and  aileron  pitching 

oscillation,  plunging  leading  the  pitching 

2 

u  -  m/-rpb  ,  mass  ratio 

oj  -  flutter  frequency 

1  /2 

u)h  -  (K^/m)  '  ,  uncoupled  plunging  natural  frequency 

1  /2 

ua  "  (Kc/I * 3a)  '  *  uncouP1ed  Patching  natural  frequency  of  airfoil 
about  elastic  axis 
1  /2 

-  (K^/I^)  >  uncoupled  pitching  natural  frequency  of  aileron 

about  elastic  axis 

(Dr  -  reference  frequency 

p  -  free  stream  air  density 

£  -  h/b,  nondimensional  plunging  degree  of  freedom 

Subscripts 

o  -  absolute  amplitude  of  displacement 

6  -  due  to  plunging  degree  of  freedom  6 

a  -  due  to  pitching  degree  of  freedom  a 

3  -  due  to  aileron  pitching  degree  of  freedom  3 
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SECTION  I 


INTRODUCTION 

The  numerical  methods  and  computer  codes  for  transonic  unsteady 
aerodynamics  have  been  developed  significantly  in  recent  years.  A 
state-of-the-art  review  of  the  numerical  methods  was  given  by  Ballhaus 
and  Bridgeman  (Reference  1).  Aeroelastic  applications  of  these 
developments  were  performed  and  reviewed  by  many  authors,  such  as 
Ashley  (Reference  2),  ftykytow  (Reference  3),  Yang,  Guruswamy,  and 
Striz  (Reference  4). 

Farmer  and  Hanson  (Reference  5}  performed  an  experimental  flutter 
analysis  of  two  dynamically  similar  wings,  one  with  supercritical 
sections  and  the  other  with  conventional  sections.  It  was  found  that 
the  experimental  results  agreed  with  those  calculated  from  the  subsonic 
lifting  surface  theory  up  to  the  Mach  number  of  0.85.  Beyond  that 
Mach  number,  the  experimental  curves  for  dynamic  pressure  showed  a 
transonic  dip  and  the  supercritical  wing  experienced  a  much  more 
pronounced  transonic  dip. 

Flutter  analyses  of  airfoils  oscillating  with  only  two  d.o.f. 's 
(plunge  and  pitch)  were  carried  out  extensively.  Rizzetta  (Reference  6) 
used  STRANS2  and  UTRANS2  (Reference  7)  to  analyze  a  NACA  64A010  airfoil. 
Yang,  Striz,  and  Guruswamy  used  STRANS2,  UTRANS2,  and  LTRAN2  (Reference  8) 
to  analyze  a  NACA  64A006  (Reference  9),  a  MBB  A-3  (Reference  10),  and 
a  NACA  64A010  (Reference  4)  airfoil.  They  also  used  STRANS2  and 
UTRANS2  to  analyze  a  TF-8A  wing  section  (Reference  4).  In  the  analysis 
of  the  streamwise  section  of  a  sweptback  wing,  Isogai  (Reference  11) 
concluded  that  the  mechanism  of  the  single  d.o.f.  flutter  dominates 
the  flutter  of  the  system  studied  at  the  bottom  of  the  transonic  dip. 
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Isogai  (Reference  12)  also  developed  a  transonic  small -perturbation 
code  which  can  be  used  for  the  reduced  frequency  kc  with  values  up 
to  1.0  and  the  entire  Mach  number  range.  He  used  the  code  to  study 
the  NACA  64A010  airfoil  for  two  cases,  simulating  characteristics 
of  a  streamwise  section  of  a  sweptback  and  an  unswept  wing,  respec¬ 
tively. 

McGrew  et^  al_.  (Reference  13)  carried  out  flutter  analysis  of 
a  TF-8A  flutter  model.  It  was  demonstrated  that  supercritical  wings 
exhibited  significantly  lower  flutter  speeds  than  a  conventional 
wing  of  equal  size  and  rigidities. 

Eastep  and  Olsen  (Reference  14)  reported  the  flutter  analysis  of 
a  rectangular  wing  by  using  the  3-D  unsteady  transonic  codes  TDSTRN 
and  TDUTRN  (Reference  7). 

In  addition  to  the  flutter  analysis,  aeroelastic  time-response 
analysis  has  also  become  a  topic  of  recent  interest.  Ballhaus  and 
Goorjian  (Reference  15)  first  performed  a  time-response  analysis  of 
a  NACA  64A006  airfoil  oscillating  with  single  pitch  d.o.f.  at  M  = 

0.88.  The  time-response  analysis  was  computed  by  using  their  program 
LTRAN2  for  unsteady  flow  coupled  with  an  integration  procedure  for 
the  structural  equation  of  motion. 

Rizzetta  (Reference  1G)  performed  a  time-response  analysis  of  a 
NACA  64A010  airfoil  with  a  single  pitch  d.o.f.  and  three  d.o.f.'s  - 
pitch,  plunge,  and  aileron  pitch.  The  LTRAN2  code  was  used.  It  was 
pointed  out  in  Reference  16  that  no  attempt  was  made  to  obtain  the 
neutrally  stable  response  curves  corresponding  to  the  flutter  condition 
for  the  three  d.o.f.  system. 

Guruswamy  and  Yang  (Reference  17)  performed  a  time-response  analysis 
of  a  NACA  64A006  airfoil  oscillating  with  plunge  and  pitch  d.o.f.'s  at 
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M  =  0.7  and  0.85,  respectively.  Again,  LTRAN2  was  used.  Parameters 
that  resulted  in  the  neutrally  stable  response  agreed  with  those 
equivalent  to  the  flutter  conditions  found  in  the  separate  flutter 
analysis.  The  print.  :ple  of  linear  superposition  of  airloads  was  used 
in  the  flutter  analysis  but  not  in  the  response  analysis. 

Recently,  several  modifications  have  been  incorporated  into 
LTRAN2.  Houwink  and  van  der  Vooren  (Reference  18)  improved  the  LTRAN2 
code  by  developing  the  NLR  (National  Aerospace  Laboratory  of  the 
Netherlands)  version.  High-frequency  terms  were  added  to  the  boundary 
conditions  and  the  pressure  computations.  Rizzetta  and  Chen  (Reference 
19)  and  Rizzetta  and  Yoshihara  (Reference  20)  included  the  cp. t  term 
in  the  smal 1 -disturbance  equation.  There  is  no  longer  the  assumption 
of  low  reduced-frequency.  Viscous  effects  were  incorporated  in  the 
code  by  using  a  viscous  ramp  method  (Reference  20). 

Goorjian  (Reference  21)  gave  a  preliminary  study  to  remove  the 
small -disturbance  and  low-frequency  restrictions  by  considering  the 
full  potential  equation.  Borland,  Rizzetta,  and  Yoshihara  (Reference  22) 
developed  a  transonic  code  LTRAN3  which  can  solve  the  3-D  low-frequency, 
unsteady  transonic  equation  by  the  time-integration  method. 

Among  these  new  developments,  the  LTRAN2-NLR  appears  to  be  an 
attractive  code  for  aeroelastic  applications.  In  the  computation  of  the 
unsteady  aerodynamic  coefficients,  the  LTRAN2-NLR  can  lift  the  limit 
of  the  reduced  frequency  (in  full  chord)  from,  say,  0.2  or  0.3  for 
LTRAN2  to  at  least  0.8.  With  this  new  capability,  the  flutter 
analysis  and  the  time-response  predictions  for  the  three  d.o.f. 

(plunge,  pitch,  and  aileron  pitch)  systems  become  more  feasible. 

More  realistic  and  broader  ranges  of  parameter  values  can  be  considered. 
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In  this  study,  LTRAN2-NLR  is  used  to  analyze  a  NACA  64A006 
conventional  and  a  MBB  A-3  supercritical  airfoil  fitted  with  trailing 
edge  ailerons  at  25%  of  the  chord.  Three  d.o.f.'s  are  considered, 
flutter  analysis  is  first  performed  and  the  effects  of  various  parameters 
are  studied.  The  neutrally  stable  time-response  predictions  are  then 
obtained.  The  effect  of  the  mass  ratio  on  the  flutter  modes  (ampli¬ 
tude  ratios  and  phase  differences)  for  various  Mach  numbers  are 
studied. 


SECTION  II 


FLUTTER  EQUATION  AND  SOLUTION  PROCEDURE 


The  parameters  and  sign  conventions  for  the  airfoil  with  pitch, 
plunge,  and  aileron  pitch  d.o.f.'s  are  defined  in  Figure  1.  The 
system  is  similar  to  that  discussed  in  Sec.  6.10  of  Reference  23. 

Based  on  the  derivations  given  in  Reference  23,  the  flutter 
equation  can  be  written  in  the  form  as, 


(1) 


p 

where  y  »  m/upb  is  the  mass  ratio;  k„  =  wc/U  is  the  reduced 

c 

frequency;  £  =  h/b  is  the  nondimensional  plunge  displacement;  the 
three  matrices  are  defined  as 
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where  r  =  (I  /mb  )  and  r0  =  ( 1,,/mb  )  are  the  radii  of  gyration  of  the 
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airfoil  and  aileron  about  the  elastic  axis  and  the  aileron  hinge 

axis,  respectively;  matrix  [A]  contains  the  9  aerodynamic  coefficients  , 

1  I 

associated  with  the  three  respective  d.o.f. 's;  =  ( K^/m) 2 ,  w  =  (K^/I^)*', 

lo„  =  (K0/I.)2  are  the  three  uncoupled  natural  frequencies;  and  u>  =  w 

p  p  p  r  a 

is  the  reference  frequency. 

The  eigenvalue  X  is  a  complex  number  defined  as 

X  =  y  (1  +  ig)  w2  b2/U2  (2) 

where  g  is  the  so-called  structural  damping  coefficient.  The  flutter 
solution  is  obtained  when  g  is  found  to  be  zero  (Reference  23).  The 
flutter  speed  is  nondimensional ized  as 

U*  =  U/boia  =  (u)r/wa)  Vm/X 

In  the  present  flutter  analysis,  the  principle  of  linear  super¬ 
position  of  airloads  was  assumed  valid  in  deriving  the  flutter  equation. 

Justification  of  this  principle  was  studied  experimentally  by 
Davis  and  Malcolm  (Reference  24)  for  a  NACA  64A010  airfoil  oscillating 
with  plunge  and  pitch  d.o.f. 's.  In  the  conclusion,  they  stated  that 
the  principles  of  superposition  and  linearity  were  shown  to  be  valid 
for  supercritical  attached  flow.  In  a  time-response  analysis  of  a 
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NACA  64A006  airfoil  oscillating  with  plunge  end  pitch  d.o.f. ' s , 

Guruswamy  and  Yang  (Reference  17)  showed  in  a  specific  example  that 
the  parameter  values  associated  with  the  flutter  speed  at  the  bottom 
of  a  transonic  dip  can  indeed  be  used  to  obtain  the  neutrally  stable 
response.  The  principle  of  superposition  was  used  in  the  flutter 
analysis  but  not  in  the  response  analysis. 

As  an  attempt  to  validate  this  principle  in  its  application  to 
flutter  analysis  of  three  d.o.f.  case,  neutrally  stable  response  results 
are  obtained  for  three  specific  examples  where  alternative  flutter 
solutions  are  found  in  the  flutter  analysis. 

In  the  flutter  analysis,  the  elastic  axis  does  not  always  coincide 
with  the  pitching  axis  for  which  the  aerodynamic  data  are  obtained. 

The  equations  for  transforming  the  aerodynamic  coefficients  for 
a  pitching  axis  O'  (see  Figure  1)  to  those  for  another  pitching  axis 
(elastic  axis)  0  are  in  the  following  form  (Reference  7) 


'16 


'la 


'16 


'la 


-  sc 


16 


CtB  C£B 


c  =  c  / 
n<5  n6 

C  a  c  1 
na  na 

cnB  =  V 


sc 


n6 


Cm6  =  W  +  S(V 


<  =  c  '  +  s(c„  ' -c  ' )  - 

\  ma  ma  'to  m6 


s2c  1 

Cf(5 


(.  cm6  =  W  +  SCm' 


(3) 


where  s  =  (ah  -  xp)/2 


SECTION  III 


EQUATIONS  OF  MOTION  AND  RESPONSE  ANALYSIS  PROCEDURE 


The  time-response  analysis  is  based  on  the  following  equations 


of  motion 


[M]  {  a" 


2  ui  2  i 

+  (mnz>  [k^“> 


r-c  i 


where  the  prime  indicates  derivative  with  respect  to  nondimensional 
time  l  -  cot. 

Equation  4  can  be  written  in  a  simplified  symbolic  form  as 


[M]  {q"}  +  [N]  {q}  =  {P} 


where  iq }  is  'e  vector  for  the  three  d.o.f. ’s;  [N]  =  (2W  /U*k  w  y 

i  c  ct 

[K] ;  and  (P >  is  the  vector  of  aerodynamic  forces.  This  equation  can 
be  solved  by  a  step-by-step  time-integration  finite  difference  approach, 
Assuming  a  linear  variation  of  acceleration,  the  velocities  and 
displacements  after  a  small  time  step  At  can  be  expressed  as 


{q,)t  *  lq,,t-A?  +  4  {q,l)t-at  +  T  (q,,)t 


<q>t  ’  +  +  4-  {q,,)t-at  +  4-  tq,,)t 


I 


Substituting  Equations  6  into  Equation  5  yields 


{q"}*  =  [F]  [ {p}^  -  [N]  (r) 


(7) 


where 

[F]  -  [CM]  [«]  (8a) 

(r)  =  <<f>t-At  +  A?{q,}E-At  +  (8b> 

This  direct  integration  method  for  structural  response  analysis 
is  well-known  (Reference  25).  In  this  study*  the  vector  for  aero¬ 
dynamic  forces  is  obtained  using  the  LTRAN2-NLR  code. 

In  each  time  step,  {q")  at  the  time  t  are  first  obtained  from  the 
known  values  of  {q},  (q ' > ,  and  {p}  at  t-At  from  Equation  6.  Based 
on  the  known  values  of  {q }  and  (q ' > .  the  aerodynamic  force  vector  {p} 
at  the  time  t  can  be  obtained  using  LTRAN2-NLR  code.  A  detailed 
description  of  the  time-response  analysis  procedure  is  given  in 
Reference  17. 

In  order  to  obtain  the  neutrally  stable  responses  for  flutter 
condition  earlier,  the  aerodynamic  equation  alone  may  be  integrated 
in  time  for  several  cycles  in  response  to  the  forced  motion.  The 
displacement  vector  of  the  airfoil  is  specified  according  to  the 
amplitude  ratios  and  the  phase  angles  of  the  flutter  mode  of  the  three 
d.o.f.  system.  After  the  aerodynamic  responses  become  periodic,  the 
system  is  set  free,  the  simultaneous  integration  procedure  begins, 
and  the  airfoil  and  the  aerodynamic  forces  drive  each  other. 
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SECTION  IV 


RESULTS  OF  AERODYNAMIC  COMPUTATIONS 

A  NACA  64A006  conventional  and  a  MBB  A-3  supe»crnic3l  airfoil 
were  studied.  Both  configurations  were  among  those  proposed  by  AGARD 
for  aeroelastic  applications  of  transonic  unsteady  aerodynamics.  The 
airfoil  coo:  Jinates  used  (Reference  26)  are  somewhat  different  from  the 
values  used  in  previous  applications.  Both  the  steady  and  unsteady  aero¬ 
dynamic  coefficients  were  computed  by  using  LTRAN2-NLR  code  (Reference  18) 

Figure  2  shows  the  steady  pressure  distributions  for  the  NACA  64A006 
airfoil  with  zero  angle  of  attack  for  M  =  0.7,  0.8,  0.85,  and  0.875.  The 
experimental  results  obtained  by  Tijdeman  (Reference  27)  at  M  =  0.85  are 
also  plotted.  At  M  =  0.85,  a  weak  shock  wave  exists  and  the  present  re¬ 
sults  agree  well  with  those  by  Tijdeman  except  in  the  neighborhood  of 
30-40%  of  the  chord.  At  M  =  0.875,  a  strong  shock  wave  is  formed. 

The  design  conditions  for  the  MBB  A-3  airfoil  are:  M  =  0.765; 
c^  =  0.58,  a  =  1.3°.  In  an  attempt  to  match  the  design  lift  coeffi¬ 
cient  of  0.58,  the  angles  of  attack  were  chosen  as  1.2°,  1.2°,  0.86°, 
and  0.6"  for  Mach  numbers  0.7,  0.74,  0.765,  and  0.78,  respectively. 

The  steady  pressure  distributions  for  tha  four  cases  were  plotted  in 
Figure  3.  Also  plotted  are  the  experimental  results  of  Bucciantini, 
Oggiano,  and  Onorato  (Reference  28)  obtained  in  Bedford  Wind  Tunnel 
for  M  =  0.765,  c  ~  0519,  and  a  =  1.5°.  For  M  =  0.765,  a  very  weak 
shock  wave  occurs  in  the  neighbornood  of  55-65%  of  the  chord.  The 
results  agree  well  with  the  experimental  data  except  in  the  neighbor¬ 
hood  of  40-60%  of  the  chord.  The  present  angle  of  attack  is,  however, 
smaller  than  that  used  in  Reference  28.  For  M  =  0.78,  the  shock  wave 
becomes  stronger  and  shifts  aft. 
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STEADY  PRESSURE  COEFFICIENT  (C 


Figure  2  Distribution  of  Steady  Pressure  Coefficients 
for  NACA  64A006  Airfoil 
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STEADY  PRESSURE  COEFFICIENT  (Cp 


-  M  -  0.700  (a  -1.2°) 

- M  -  0.740  (a  *  1 . 2° ) 

-  M  =  0.765  (  a  -0.86°) 

.  M  -  0  7A£71  (  a  *01  R0°1 


X  /  C 


Figure  3  Distribution  of  Steady  Pressure  Coefficients 
for  MBB  A-3  Airfoil  with  coai0.58 


In  the  unsteady  aerodynamic  calculations,  60  time  steps  were  used 
for  each  half-cycle.  The  amplitudes  chosen  were  0.1C  and  0.1°  for 
the  plunge  and  pitch  d.o.f.'s,  respectively.  The  nine  aerodynamic 
coefficients  cAfi.  c^,  cm(S,  c^.  cmg»  cn6,  cnQt.  and  cng  were 
computed  by  using  a  79  x  99  grid.  Both  airfoils  fitted  with  trailing 
edge  ailerons  of  25%  chord,  were  assumed  to  pitch  about  the  quarter 
chord  axis.  The  results  were  given  in  Tables  1  to  8  with  reduced 

frequency  k  equal  to  0.1,  0.2,  0.3,  0.4,  0.6,  0.8,  1.2,  and  1.6  for 

v 

each  Mach  number  considered. 
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Table  1.  Aerodynamic  Coefficients  for  the  NACA  64A006 


Airfoil  Pitching  at  the  Quarter  Chord  at  M=0.7 


1 

Cjz  s 

Cj2« 

Cj>a 

CmS 

^mo 

1 

Real 

I  mag 

Real 

Imag 

Rea  1 

Imag 

Real 

Imag 

Real 

Imag 

.05 

.066 

0.41 

8.28 

-1  .31 

5.14 

-0.82 

.005 

-.010 

-.205 

-.091 

0.1 

.161 

0.76 

7.55 

-1  .61 

4.62 

-1 .24 

.015 

-.019 

-.189 

-.153 

0.2 

.385 

1  .24 

6.44 

-1  .73 

3.92 

-1 .28 

.028 

-.042 

-.166 

-.287 

0.3 

.508 

1  .65 

5.91 

-.940 

3.43 

-1 .32 

.061 

- .  063 

-.160 

-.416 

0.4 

.557 

2.02 

5.54 

-.290 

3.14 

-1  .21 

.108 

-.085 

-.148 

-.553 

0.5 

.520 

2.46 

5.28 

0.01  2 

2.93 

-1.12 

.176 

-.108 

-.146 

-.686 

0.6 

.491 

2.73 

5.12 

0.538 

2.77 

-1  .06 

.245 

-.135 

-.130 

-.822 

0.8 

.255 

3.47 

4.95 

1  .  327 

2.60 

-0.84 

.427 

-.203 

-.058 

-1.11 

1  .2 

-.77 

5.03 

4.92 

3.140 

2.41 

-0.60 

1  .03 

-.365 

-.110 

-1  .69 

1  .6 

-1  .9 

7.06 

5.27 

4.530 

2.24 

-0.52 

1  .79 

-.715 

-.178 

-2.37 

| 

Cm  B 

Cn  f 

cr 

(X 

Cne 

fl 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

.05 

m 

.04 

-0.027 

0.0001 

-.0025 

-.0502 

- . 0005 

-.0879 

-.0018 

0.1 

.03 

-0.036 

0.0002 

- . 0047 

-.0467 

-.0025 

-.0856 

-.0019 

0.2 

.03 

-0.063 

0 . 0005 

-.0084 

-.0416 

-.0135 

-.0831 

-.0059 

0.3 

-i 

.03 

-0.088 

0.0023 

-.0118 

-.0376 

-.0244 

-.0813 

-.0114 

0.4 

-i 

.04 

-0.114 

0.0053 

-.0151 

-.0346 

-.0346 

-.0806 

-.0171 

-i 

.06 

-0.137 

0.0093 

-.0187 

-.0323 

- . 0445 

-.0805 

-.0218 

0.6 

-i 

.07 

-0.160 

0.0148 

-.0215 

-.0312 

-.0539 

-.0806 

- .0272 

0.8 

-i 

.09 

-0.196 

0.0280 

-.0280 

-.0205 

-.0754 

-.0809 

-.0364 

1 .2 

-i 

.16 

-0.252 

0.0720 

-.0400 

-.0045 

-.1210 

-.0829 

-.0568 

1 .6 

-i 

.21 

-0.269 

0 . 1  350 

-.0550 

0.0146 

-.1710 

-.0860 

-.0748 
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Table  2.  Aerodynamic  Coefficients  for  the  NACA  64A006 

Airfoil  Pitching  at  the  Quarter  Chord  at  M=0.8 


Cj2  s 

c, 

2a 

Cj Za 

Cm  s 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

.05 

.106 

0 

50 

1C 

).0 

-2.13 

6 

.12 

-1 

.  64 

.006 

—  t 

01  3 

— 

253 

—  # 

129 

0 

.  1 

.278 

0 

86 

8. 

57 

-2.78 

5 

.20 

-1 

.99 

.020 

- 

025 

- 

252 

- . 

204 

0 

.2 

.565 

1 

33 

7. 

02 

-2.28 

4 

.04 

-2 

.06 

.037 

- 

061 

- 

251 

- . 

386 

0 

.3 

.750 

1 

70 

6. 

23 

-1  .67 

3 

.  39 

-1 

.96 

.075 

- 

095 

- 

281 

552 

0 

.4 

.818 

2 

08 

5. 

84 

-.925 

3 

.02 

-1 

.75 

.143 

- 

134 

- 

31  8 

- . 

713 

0 

.5 

.810 

2 

42 

5. 

57 

-.585 

2 

.75 

-1 

.59 

.210 

- 

190 

- 

344 

- . 

896 

0 

.6 

.803 

2 

76 

5. 

41 

-.094 

2 

.48 

-1 

.55 

.312 

- 

256 

- 

416 

-1 

.09 

c 

.8 

.707 

3 

50 

5. 

36 

0.563 

2 

.10 

-1 

.37 

.523 

- 

419 

- 

548 

-1 

.43 

1 

.2 

.417 

4 

98 

5. 

10 

1  .590 

1 

.47 

-0 

.85 

.857 

- 

952 

- 

968 

-1 

.79 

1 

.6 

-.45 

6 

28 

4. 

76 

2.560 

1 

.39 

-0 

.26 

1.18 

-1 

.36 

-1 

.07 

-1 

.91 

CmS 

cr 

s 

cr 

cx 

Cna 

kc 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Rea  1 

Imag 

.05 

-1 

28 

-0.019 

0.0000 

-.0025 

-.0505 

-.0001 

-.0960 

- . 0004 

0.1 

-1 

28 

-0.032 

0.0002 

-.0046 

-.0462 

-.0024 

-.0933 

-.0023 

0.2 

-1 

30 

-0.044 

0.0010 

-.0084 

-.0415 

-.0159 

-.0915 

-.0085 

0.3 

-1 

33 

-0.047 

0 . 0045 

-.0115 

-.0376 

-.0305 

-.0911 

-.0152 

0.4 

-1 

36 

-0.039 

0.0078 

-.0155 

-.0347 

-.0429 

-.0920 

-.0220 

0.5 

-1 

.39 

-0.019 

0.0130 

-.0190 

-.0352 

-.0542 

-.0935 

-.0270 

0.6 

-1 

41 

0.0170 

0.0193 

-.0239 

-.0306 

-.0688 

-.0950 

-.031  8 

0.8 

-1 

.42 

0 . 1  1  60 

0.0369 

-.0325 

-.0259 

-.0966 

-.0976 

-.0415 

1 .2 

-1 

19 

0.2910 

0.0842 

-.0572 

-.0237 

-.1510 

-.1026 

-.0539 

1 .6 

-0 

.93 

0 . 1  450 

0 . 1 450 

-.0872 

-.0165 

-.1.10 

-.1014 

-.0688 
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Table  3.  Aerodynamic  Coefficients  for  the  NACA  64A006 


Airfoil  Pitching  at  the  Quarter  Chord  at  M=0.85 


s 

Cj 

2  ot 

Cj2e 

S 

Cmot 

1 

8 

Real 

I  mag 

Real 

I  mag 

Rea  1 

Imag 

Real 

Imag 

Real 

Imag 

.05 

.210 

0 

61 

i: 

1.2 

-4.20 

7 

.41 

-2 

85 

.006 

—  ( 

026 

-.510 

-.118 

0 

.  1 

.432 

0 

97 

9. 

72 

-4.33 

5 

.65 

-3 

26 

.026 

- 

051 

-.505 

-.257 

0 

.2 

.843 

1 

37 

7. 

40 

-3.45 

3 

.85 

-3 

00 

.032 

- 

1  29 

-.584 

-.482 

0 

.3 

1  .05 

1 

65 

6. 

34 

-2.69 

2 

.85 

-2 

66 

.053 

- 

210 

-.714 

-.666 

0 

.4 

1  .22 

1 

90 

5. 

86 

-2.13 

2 

.28 

-2 

28 

.085 

- 

316 

-.888 

-.772 

0 

.5 

1  .22 

2 

18 

5. 

32 

-1  .03 

1 

.82 

-1 

82 

.081 

- 

382 

-.995 

-.778 

0 

.6 

1  .22 

2 

39 

4. 

84 

-1.12 

1 

.66 

-1 

34 

.073 

- 

459 

-1  .01 

-.707 

0 

.8 

.974 

2 

96 

4. 

48 

0.003 

1 

.68 

-0 

86 

.165 

- 

517 

-.876 

-.846 

1 

.2 

.305 

4 

41 

4. 

57 

1  .364 

1 

.72 

-0 

65 

.570 

- 

790 

-.366 

-1  .31 

1 

.6 

-.49 

5 

92 

4. 

58 

2.865 

1 

.58 

-0 

43 

1  .04 

-1 

.22 

-.939 

-1  .72 

| 

CmB 

cr 

s 

cr 

c X 

Cn  6 

13 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

in 

o 

-1  .66 

0.0300 

0.0002 

-.0019 

- . 0375 

-.0039 

-.  1000 

-.0035 

0.1 

-1  .67 

0.0300 

0.0012 

- . 0035 

-.0351 

-.0121 

-.0938 

-.0082 

0.2 

-1  .72 

0.0900 

0.0032 

-.0072 

-.0340 

-.0275 

-.0944 

-.01 84 

SSkj 

-1  .76 

0 . 2470 

0.0075 

-.0115 

-.0332 

-.0425 

-.0997 

-.0249 

0.4 

-1  .70 

0 . 4560 

0.0135 

-.0165 

-.0341 

-.0590 

-.  1060 

-.0304 

0.5 

-1  .51 

0.5810 

0.0190 

-.0230 

-.0360 

-.0740 

-.1110 

-.0317 

0.6 

-1  .26 

0 . 5870 

0.0245 

-.0306 

-.0410 

-.0840 

-.1120 

-.0300 

0.8 

-1  .04 

0 . 3450 

0.0386 

-.0430 

-.0400 

-.1040 

-.1080 

-.0350 

1  .2 

-.974 

0.1660 

0.0845 

-.0670 

-.0344 

-.1541 

-.1100 

-.0520 

1  .6 

-.961 

0.0974 

0.1  460 

-.1020 

-.0230 

-.2060 

-.1130 

-.0670 

Table  4.  Aerodynamic  Coefficients  for  the  NACA  64A006 


Airfoil  Pitching  at  the  Quarter  Chord  at  M=0.875 


1 

1 

Cj 

2£ 

— 

Cj 

26 

Cn 

nS 

Cmot 

1 

1 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

.05 

.666 

0 

69 

14 

.0 

-12.7 

7.91 

-8 

37 

-.15 

—  ^ 

143 

-2.89 

2. 

733 

0. 

1 

.954 

0 

81 

8. 

25 

-8.72 

4.41 

-6 

.07 

-.21 

- . 

157 

-1  .64 

1  . 

874 

0  . 

2 

1.14 

1 

03 

5. 

78 

-5.33 

2.36 

-3 

.64 

-.26 

- . 

167 

-1  .03 

1  . 

089 

0. 

3 

1.18 

1 

27 

4. 

93 

-3.55 

1  .83 

-2 

.39 

-.26 

169 

-.842 

0. 

619 

0. 

4 

1.16 

1 

53 

4. 

47 

-2.35 

1  .75 

-1 

.64 

-.23 

169 

-.699 

0. 

288 

0. 

5 

1.11 

1 

86 

4. 

36 

-1  .55 

1  .75 

-1 

.29 

-.16 

191 

-.607 

0. 

014 

0. 

6 

1  .07 

2 

21 

4. 

36 

-0.55 

1  .81 

-1 

.09 

-.08 

- . 

243 

-.605 

- . 

252 

0. 

8 

.880 

2 

87 

4. 

40 

0.057 

1  .75 

-0 

.92 

.113 

367 

-.688 

- . 

653 

1 . 

2 

.350 

4 

18 

4. 

25 

1  .203 

1  .58 

-0 

.54 

.490 

736 

-.848 

-1 

.15 

1 

6 

-.40 

5 

.61 

4. 

27 

2. 177 

1  .52 

-0 

.34 

.925 

-1 

.15 

-.905 

-1 

.54 

■ 

CmB 

cr 
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Table  7.  Aerodynamic  Coefficients  for  the  MBB  A-3  Airfoil 


at  M=0.765  and  a^O.So0 
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Table  8.  Aerodynamic  Coefficients  for  the  MBB  A-3  Airfoil 


at  M=0.78  and  a=0.6° 


■ 

wm 

Ci 

2a 

C  J 2e 

S 

Cm  cr 

m3 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

0.1 

.634 

.981 

10.2 

-5.63 

5.76 

-3.86 

-.07 

-.117 

-1  .22 

.4544 

0.2 

1  .01 

1  .31 

7.32 

-4.34 

3.74 

-3.22 

-.12 

-.163 

-.917 

.2452 

0.3 

1 .17 

1  .55 

6.06 

-3.14 

2.87 

-2.52 

-.13 

-.1  83 

-.783 

.0354 

0.4 

1  .20 

1  .80 

5.46 

-2.21 

2.49 

-2.01 

-.11 

-.192 

-.707 

-.150 

0.6 

1  .06 

2.39 

4.91 

-.773 

2.19 

-1 .34 

.009 

-.229 

-.614 

-.498 

0.8 

.809 

3.06 

4.73 

0.287 

2.08 

-1  .01 

.210 

-.315 

-.584 

-.832 

1  .2 

.219 

4.62 

4.83 

1  .743 

2.00 

-.736 

.690 

-.638 

-.643 

-1  .46 

1  .6 

-.62 

6.26 

5.02 

2.717 

1  .82 

-.574 

1  .22 

-1.15 

-.811 

-1  .97 

| 

Cm/3 

C„s 

cr 

cx 

C  n  /3 

B 

Real 

Imag 

Rea  1 

Imag 

Rea  1 

Imag 

Rea  1 

Imag 

0.1 

-1  .85 

0.4989 

.00491 

-.0002 

.00378 

-.0487 

-.0584 

-.0342 

0.2 

-1  .57 

0.5260 

.01081 

- . 0052 

-.0166 

- . 0595 

-.0766 

-.0429 

0.3 

-1  .38 

0.4967 

.01540 

-.0120 

-.0279 

-.0664 

-.0886 

-.0401 

0.4 

-1  .25 

0.4378 

.01939 

-.0194 

-.0337 

-.0731 

-.0947 

-.0385 

0.6 

-1  .09 

0.2808 

.02898 

-.0324 

-.0377 

-.0869 

-.0979 

-.0363 

0.8 

-1  .05 

0.1559 

. 04245 

-.0437 

-.0361 

-.1045 

-.0978 

-.0408 

1 .2 

-1  .08 

0.0156 

.08425 

-.0653 

-.0232 

-.1468 

-.0951 

-.0568 

1  .6 

-1  .08 

-.0036 

.14710 

-.0962 

-.0092 

_ 

-.1964 

-.0967 

-.0734 

22 


SECTION  V 


RESULTS  OF  FLUTTER  ANALYSIS 

Unless  otherwise  specified,  the  present  flutter  analysis  was 
based  on  the  following  set  of  common  parameters:  a^  =  -0.2;  xa  =  0.2; 
r  =  0.5;  xQ  =  0.008;  r0  =  0.06;  w./u>  =  0.3;  and  y  =  50.  The  Mach 
numbers  considered  were  0.7,  0.8,  0.85,  and  0.875  for  the  NACA  64A006 
airfoil  and  0.7,  0.74,  0.765,  and  0.78  for  the  MBB  A-3  airfoil.  A 
flutter  analysis  of  a  flat  plate  with  a  trail ing-edge  aileron  at  M  =  0.7 
was  first  performed  using  the  aerodynamic  data  obtained  from  the  kernel 

J 

function  method  and  LTRAN2-NLR  code.  The  two  sets  of  curves  obtained 
for  U*  vs.  by  the  two  methods  agreed  well  with  each  other 

(Reference  29). 

1.  Effect  of  Frequency  Ratios  (wg/ooa  and  u^/^) 

The  curves  for  the  flutter  speed  vs.  %/a)  ranging  from  0  to  2 

p  a 

for  three  values  (0.1,  0.3,  and  0.5)  and  various  Mach  numbers 

were  obtained  for  both  airfoils  and  plotted  in  Reference  29.  Only  the 
results  for  the  NACA  64A006  airfoil  at  M  =  0.85  is  given  here  (Figure  41. 

The  flutter  boundaries  for  the  three  branches:  bending-torsion; 
bending-aileron;  and  torsion-aileron  can  be  correlated  to  those  for  the 
corresponding  binary  cases  except  in  the  transition  zones.  Such 
correlation  was  made  in  the  study  of  subsonic  case  (Reference  30). 

For  w^/w^  =  0.1  in  Figure  4,  the  flutter  speed  for  the  bending- 

torsion  branch  increases  very  slightly  with  u)Q/w  and  the  curve  approaches 

p  a 

asymptote  for  w  /to  -*»,  which  corresponds  to  the  bending-torsion  binary 

p  a 

flutter  value.  For  w  /w  =  o.l,  the  flutter  boundaries  for  the  torsion- 

h  a 

aileron  and  bending-aileron  branches  appear  in  the  regions  of  w <  0.92 
and  0.18,  respectively. 


23 


REDUCED  FREQ.  (Kc)  FLUTTER  SPEED  (U/bU^I 


NRCR  64R006  RIRFOIL  RT  M-0.85 


RILERON-RIRFOIL  PITCH  FREQUENCY  RRTIO  (COp/OJ*) 


Figure  4  Effect  of  Aileron-Airfoil  Pitch  Frequency  on 
Flutter  Speed  for  Three  Different  Airfoil 
Plunge -Pitch  Frequency  Ratios 


For  oju/(o  =  0.3  and  0.5,  the  flutter  boundaries  for  the  torsion- 
h  a 

aileron  branches  are  not  too  much  different  than  those  for  u)^/wa  =  0.1. 
However,  the  two  boundaries  for  bending-torsion  and  bending-aileron 
branches  become  continuous  curves,  each  with  a  transition  portion. 

The  resulting  curves  obtained  for  all  other  Mach  numbers  for  both  air¬ 
foils  show  the  trends  similar  to  those  found  in  Figure  4  of  Reference  29. 

A  common  means  to  eliminate  the  flutter  boundaries  for  the 
bending-aileron  and  the  torsion-aileron  branches  is  to  use  aileron 
mass  balance  or  high  aileron  pitch  stiffness.  For  a  modern  aircraft 
with  power-operated  irreversible  controls,  the  value  of  is 

much  higher  than,  say,  1.0  and  thus  the  flutter  boundaries  for  the 
two  aileron  associated  branches  do  not  exist.  In  the  following  analyses, 

a)  /t jo  values  considered  were  0.8  and  1.5,  respectively,  such  that  the 

0  a 

flutter  boundaries  for  the  bending-aileron  branch  is  totally  avoided 
but  the  torsion-aileron  branch  still  exists  atwg/u)a  =  0.8. 

2.  Effect  of  Mass  Ratio  (y) 

The  curves  for  the  flutter  speed  vs.  mass  ratio  y  for  two  u)g/o>a 

values,  0.8  and  1.5,  and  various  Mach  numbers  were  obtained  for  both 

airfoils  (Reference  29).  Only  the  results  for  the  MBB  A-3  airfoil 

at  the  design  Mach  number  0.765  is  shown  here  (Figure  5). 

Foru^/w  =  1.5,  the  flutter  speed  for  the  bending-torsion  branch 
p  a 

increases  with  y  in  a  fashion  similar  to  that  for  the  bending-torsion 
binary  flutter.  The  flutter  boundaries  for  the  two  aileron  associated 
branches  are  absent. 
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MBB-A3  AIRFOIL  AT  M-0.765 


For  u)„/<jj  =  0.8,  the  flutter  speed  for  the  bending-torsion 
p  a 

branch  increases  with  n  at  a  much  slower  rate  than  that  in  the 
former  case,  due  to  the  effect  of  aileron  motion.  The  boundary 
for  the  bending-torsion  branch  is  lower  than  the  torsion-aileron 
branch. 

3.  Effect  of  the  Position  of  Airfoil  Mass  Center  (x^) 

The  curves  for  the  flutter  speed  vs.  ranging  from  -0.2  to 
0.4  for  two  u)0/w  values  (0.8  and  1.5}  were  obtained  for  the  MBB 
A-3  and  the  NACA  64A006  airfoils  at  M  =  0.765  and  0.85,  respectively 
(Reference  29).  Only  the  results  for  the  MBB  A-3  airfoil  at  the 
design  Mach  number  0.765  is  shown  here  (Figure  6). 

For  =  1.5,  the  flutter  speed  for  the  bending-torsion  branch 

increases  as  the  mass  center  moves  forward.  The  flutter  boundary 
for  the  torsion-aileron  branch  only  appears  when  xa  <  0.  The  flutter 
boundary  for  the  bending-aileron  branch  is  absent. 

For  oj0/ oj  =  0.8,  the  flutter  boundary  for  the  torsion-aileron 

p  ex 

branch  extends  to  x^  =  0.12.  The  boundary  for  the  bending-torsion 
branch  appears  at  x^  >  -0.1.  The  last  two  boundaries  become  continuous 
with  a  transition  portion  in  the  neighborhood  of  xa  =  0.1. 

It  is  physically  clear  that  when  an  airfoil  oscillates  with  only 
a  single  pitching  motion,  the  forward  movement  of  the  mass  center 
stabilizes  the  airfoil.  Such  stabilizing  effect  is  also  evident  in 
the  present  flutter  speed  curves  for  the  bending-torsion  branch  of  a 
three  d.o.f.  system. 
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4.  Effect  of  the  Position  of  Elastic  Axis  (a.  ) 

h 

The  curves  for  the  flutter  speed  vs.  ah  ranging  from  -0.5  to 

0  for  two  (*)./<*)  values  (0.8  and  1.5)  were  obtained  for  the  NACA 
p  a 

64A006  and  the  MBB  A-3  airfoils  at  M  =  0.85  and  0.765,  respectively 
(Reference  29).  Only  the  results  for  the  NACA  64A006  airfoil  is 
shown  here  (Figure  7). 

For  ajg/co^  =  1.5,  the  flutter  speed  for  the  bending-torsion  branch 
increases  as  the  elastic  axis  moves  forward.  The  flutter  boundary 
for  the  torsion-aileron  branch  appears  in  the  form  of  a  loop  in  the 
region  of  ah  <  -0.47. 

Foroig/co^  =  0.8,  the  flutter  boundary  for  the  bending-torsion 
branch  assumes  a  similar  trend  as  that  of  the  former  case.  The  loop¬ 
shaped  boundary  for  the  torsion-aileron  branch  extends  to  ah  =  -0.14. 
The  boundaries  for  the  bending-aileron  branch  disappear  in  both  cases. 

The  forward  movement  of  the  elastic  axis  has  the  same  stabilizing 
effect  as  that  of  the  mass  center.  Thus,  the  trend  of  the  flutter 
boundaries  of  the  bending-torsion  branch  observed  in  Figure  7  is  quite 
similar  to  that  in  Figure  6. 

5.  Effect  of  the  Position  of  Aileron  Mass  Center  (xD) 

p 

The  curves  for  the  flutter  speed  vs.  xQ  ranging  from  0  to  0.02 

p 

for  two  o)g/wa  values  (0.8  and  1.5)  were  obtained  for  the  NACA  64A006 
and  the  MBB  A-3  airfoils  at  M  =  0.85  and  0.765,  respectively  (Reference 
29).  Only  the  results  for  the  NACA  64A006  airfoil  are  shown  here 
(Figure  8). 

For  w0/(jo  =  1.5,  only  the  flutter  boundaries  for  the  bending- 
p  a 

torsion  branch  appear.  The  flutter  speed  decreases  slightly  as 
xQ  increases. 

P 


.-■s?  n't'  mVlfr'ir. 
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REDUCED  FREQ.  <Kcf  FLUTTER  SPEED  (U/bOOJ 


Figure  7  Effect  of  Position  of  Elastic  Axis  on  Flutter  Speed 
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REDUCED  FREQ.  (Kc)  FLUTTER  SPEED  (U/bUy 


For  =  0.8,  the  decrease  in  flutter  speed  for  the  bendinq- 

torsion  branch  with  x(i  is  more  pronounced  than  that  of  he  former  case. 
The  flutter  boundary  for  the  torsion-aileron  branch  appears  only  at 
Xg  >  0.004.  The  flutter  boundary  for  the  bending-aileron  branch  is 
absent. 

It  may  be  seen  that  the  forward  movement  of  the  aileron  mass 
center  can  not  only  eliminate  the  flutter  boundary  for  the  torsion- 
aileron  branch  at  x0  <  0.0025  in  this  case,  but  also  increase  the  flutter 

p 

speed  of  the  bending-torsion  branch.  The  advantage  of  using  aileron 
mass  balance  is  seen  in  this  case. 

6.  Effect  of  Mach  Number  (M) 

Results  for  the  flutter  speed  vs.  Mach  number  were  plotted  for 
both  airfoils  for  y  =  20  and  50  in  Figures  9.  The  curves  were  plotted 

for  two  values  of  u)0/io  :  0.8  and  1.5.  All  these  curves  are  for  the 

p  a 

bending-torsion  branch  only. 

For  oj0/ u>  =1.5  and  for  the  present  values  assumed  for  the  other 
aeroelastic  parameters,  only  the  flutter  boundaries  for  the  bending- 
torsion  branch  are  present.  Those  for  the  other  two  aileron  related 

branches  are  absent.  For  w./w  =0.8,  the  flutter  boundaries  for  the 

P  a 

bending-aileron  branch  are  absent  and  those  for  the  torsion-aileron 

branch  are  lower  than  those  for  the  bending-torsion  branch.  However, 

the  flutter  speed  curves  for  the  bending-torsion  branch  at  u>0/w  =  0.8 

p  ot 

were  still  plotted  in  Figure  9  for  the  purpose  of  comparison  with  those 

at  u>  / u)  =1.5.  As  noted  earlier,  for  a  modern  aircraft  with  power- 
p  a 

operated  irreversible  controls  the  value  of  w  /u>a  is  much  higher  than 

1.0. 
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FLUTTER  SPEED  (U/bCJa)  FLUTTER  SPEED  (U/bCOa) 


0.6  0.65  0.7  0.75  0.8  0.85  0.8 


MACH  NUMBER  (  M  ) 


Figure  9  Effect  of  Mach  Number  on  Flutter  Speed  for 
Both  Airfoils  for  Two  Mass  Ratios: 

(a)  y=20  and  (b)  y=50 
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The  figure  shows  that  the  bottom  of  the  transonic  dip  occurs 
in  the  neighborhood  of  the  design  Mach  number  of  0.765  for  the  MBB  A-3 
airfoil  and  in  the  neighborhood  of  M  =  0.85  for  the  NACA  64A006  airfoil. 
In  a  flutter  analysis  of  the  NACA  64A006  airfoil  with  only  plunge  and 
pitch  d.o.f. 's  (Reference  9),  the  transonic  dip  for  the  parameter 
values  considered  was  also  found  to  be  near  M  -  0.85. 

The  figure  shows  that  the  mass  ratio  has  a  detrimental  effect 
in  deepening  the  transonic  dip.  This  effect  was  pointed  out  in  Reference 
3. 

The  present  examples  show  that  the  transonic  dips  for  the  MBB 
A-3  supercritical  airfoil  are  considerably  deeper  and  occur  at  earlier 
Mach  numbers  than  thc^e  for  the  NACA  64A006  airfoil.  It  should,  however, 
be  noted  that  both  airfoils  are  not  comparable  in  configuration. 

Results  for  flutter  speed  vs.  Mach  number  for  other  values  of  the 
parameters  (wh/u>a>  x^,  a^)  are  available  in  Reference  29. 

7.  Discussion  of  the  Flutter  Modes 

To  study  the  flutter  mode,  the  case  of  two  d.o.f. 's  (plunge  and 
pitch)  was  first  considered.  The  flutter  speed,  plunge-pitch  amplitude 
ratio,  and  plunge-pitch  phase  difference  were  plotted  against  the  mass 
ratio  in  Figure  10  for  both  airfoils  and  various  Mach  numbers.  The 
values  of  the  parameters  used  are  the  same  as  those  assumed  for  a 

NACA  64A010  airfoil  in  Case  A  of  Reference  12:  x  =  1.8;  r  =  1.867; 

a  a 

a,  =  -2.0;  and  w. /w  =  1.0. 
h  h  a 

It  is  seen  that,  at  the  bottoms  of  the  transonic  dips  (M  1  0.78 
and  0.875  for  MBB  A-3  and  NACA  64A006  airfoils,  respectively),  the 
amplitude  ratios  gradually  converge  to  a  constant  value  of  1.868 
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(first  natural  mode)  and  the  phase  differences  gradually  approach 
zero  as  the  mass  ratio  becomes  larger.  The  phenomenon  observed  here 
agrees  with  that  pointed  out  in  Reference  12. 

As  shown  in  Figure  11,  for  the  NACA  64A006  airfoil,  the  two  flutter 
speed  curves  for  M  =  0.85  and  0.875  cross  each  other  at  y  =  185.  For 
higher  y  values,  the  flutter  speed  for  M  =  0.875  becomes  increasingly 
lower  than  that  for  M  =  0.85  and  the  bottom  of  the  transonic  dip 
definitely  shifts  to  a  Mach  number  higher  than  0.85. 

It  is  seen  that  for  M  =  0.875  and  y  >  185,  the  amplitude  ratio 
approaches  a  constant  value  of  12.86  (first  natural  mode)  and  the  phase 
difference  decreases  to  zero  as  y  becomes  larger. 

For  the  MBB  A-3  airfoil,  the  two  flutter  curves  for  M  =  0.765 
and  0.78  cross  each  other  at  y  =  195.  For  M  =  0.78  and  y  >  195,  the 
amplitude  ratio  approaches  rather  slowly  the  constant  value  of  12.86 
(first  natural  mode)  and  the  phase  difference  decreases  to  zero  as  y 
becomes  larger. 

The  MBB  A-3  airfoil  with  three  d.o.f.’s  is  then  considered.  The 
parameter  values  are  defined  in  the  first  paragraph  of  this  section. 

It  is  also  defined  thatu0Aj  =  1.5.  The  results  are  given  in  Figure 

pa 

12. 

At  M  =  0.78,  the  plunge-pitch  and  the  plunge-aileron  pitch 
amplitude  ratios  approach  two  constant  values  of  12.81  and  10.1  respectively 
(first  natural  mode)  and  the  two  phase  differences  approach  zero  as 
y  becomes  larger. 
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SECTION  VI 


RESULTS  OF  TIME  RESPONSE  ANALYSIS 

Time  response  results  were  obtained  for  the  NACA  64A006  airfoil 
at  M  =  0.7  and  0.85  and  for  the  MBB  A-3  airfoil  at  M  =  0.765, 
respectively.  The  parameters  are  defined  as:  ah  =  -0.2,  =  0.2; 

r  =  0.5;  x.  =  0.008;  r 0  =  0.06;  uj./cj  =  0.3;  w0/m  =  1.5;  and  k  =  0.3. 

a  B  8  h  a  B  a  c 

A  time-response  analysis  was  first  performed  for  the  NACA 
64A006  airfoil  at  M  =  0.7.  Based  on  the  flutter  speed  found  in  a  flutter 
analysis,  U*  =  2.816  at  y  =  48,  the  neutrally  stable  responses  were 
indeed  obtained  (Reference  29). 

A  time-response  analysis  was  then  performed  for  the  NACA  64A006 
airfoil  at  M  =  0.85.  The  results  for  the  three  displacements  and  the 
three  aerodynamic  forces  were  plotted  in  Figures  13  und  14,  respectively. 
In  Figure  14,  c^,  cmQ,  and  cnQ  are  the  amplitudes  of  c^, 

cm*  and  cn*  resPectively,  obtained  in  the  forced  harmonic  motion. 

The  airfoil  was  first  forced  to  oscillate  for  5-1/2  cycles  ip 
order  for  the  response  of  the  three  aerodynamic  coefficients  to  become 
periodic.  The  airfoil  was  then  set  free  and  the  aeroelastic  time- 
responses  were  calculated.  Based  on  a  flutter  solution,  U*  =  2.436 
at  y  =  52.1,  the  time-responses  were  found  to  be  slightly  diverging 
(unstable).  The  flight  speed  was  then  reduced  by  10%  and  the  time- 
responses  were  found  to  be  converging  (stable).  When  the  flight 
speed  was  reduced  by  only  4%,  the  neutrally  stable  responses  were 
eventually  obtained.  In  the  neutrally  stable  free  motion,  the  frequency 
was  found  to  be  1%  higher  than  that  of  the  original  forcing  function. 

For  the  MBB  A-3  airfoil  at  M  =  0.765,  the  response  results  for 

the  three  displacements  and  the  three  aerodynamic  forces  were  plotted 

* 

in  Figures  15  and  16,  respectively.  In  Figure  16,  dc^  is  the  differential 
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Figure  14  Effect  of  Flutter  Speed  on  Aerodynamic  Responses 
for  the  NACA  64A006  Airfoil  at  M  =  0.85 
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relative  to  the  mean  value  of  c^  in  forced  motion;  dc^  and 

dcn  are  defined  in  a  similar  way. 

Again,  the  airfoil  was  first  forced  to  oscillate  for  5  and  1/2 
cycles  to  obtain  the  periodic  aerodynamic  responses.  Based  on  a 
flutter  solution,  U*  =  2.456  at  y  =  58,  the  time-responses  were  found 
to  be  slightly  diverging.  The  neutrally  stable  responses  were  obtained 
when  the  flight  speed  was  reduced  by  3.5%. 


SECTION  VII 


CONCLUDING  REMARKS 

Based  on  the  present  flutter  and  time-response  analyses,  the 
following  concluding  remarks  can  be  made. 

(1)  As  can  be  seen  in  Figures  2  and  3,  strong  shocks  are 
present  at  M  s  0.875  and  0..78  for  the  NACA  64A006  and  MBB  A-3  airfoils, 
respectively.  Both  the  steady  and  unsteady  aerodynamic  results  were 
obtained  by  using  the  LTRAN2-NLR  code.  Difficulty  was,  however, 
encountered  when  attempting  to  obtain  unsteady  results  for  the  NACA 
64A006  airfoil  at  M  =  0.9. 

(2)  In  computing  the  aerodynamic  coefficients  for  all  cases, 
k  values  of  up  to  1.6  were  considered.  For  a  flat  plate  at  M  =  0.7, 

w 

all  the  nine  unsteady  aerodynamic  coefficients  were  obtained  which 
agreed  reasonably  well  with  the  kernal  function  results  (Reference  29). 

(3)  For  the  parameter  values  considered,  the  results  show  that 
the  flutter  boundaries  for  the  bending-aileron  and  the  torsion-aileron 
branches  generally  disappear  when  w0/u>  is  greater  than  one.  The  aileron 

P  ct 

effect  on  flutter  speed  of  the  bending- torsion  branch  diminishes 
gradually  as  increases. 

(4)  As  seen  in  Figures  5,  6,  and  7,  the  flutter  speed  for 

w  =  1.5  increases  as  the  mass  ratio  becomes  larger,  as  the  mass 
center  moves  forward,  and  as  the  elastic  axis  moves  forward,  respectively. 
It  also  increases  as  the  radius  of  gyration  of  the  airfoil  becomes 
larger  (Reference  29). 
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(5)  As  seen  in  Figure  8,  forward  movement  of  the  aileron  mass 
center  can  not  only  eliminate  the  flutter  boundaries  of  aileron 
associated  branches  but  also  increase  the  flutter  speed  of  the 
bending-torsion  branch. 

(6)  For  the  sets  of  parameter  values  considered  and  for  the 
NACA  64A006  and  MBB  A-3  airfoils  at  M  =  0.875  and  0.78,  respectively, 
each  flutter  mode  converges  to  its  corresponding  first  natural  mode 
as  the  mass  ratio  becomes  larger. 

(7)  For  the  NACA  64A006  airfoil  at  M  =  0.7,  neutrally  stable 
responses  were  obtained  based  on  a  set  of  parameter  values  correspond¬ 
ing  to  a  flutter  condition  in  a  separate  flutter  analysis  (Reference 
29).  For  the  NACA  64A006  airfoil  at  M  =  0.85  and  the  MBB  A-3  airfoil 
at  M  =  0.765,  neutrally  stable  responses  were  obtained  based  on  the 
flight  speed  4%  and  3.5%  lower  than  the  respective  flutter  speeds 
found  in  the  separate  flutter  analyses.  Such  small  differences  may 
mainly  be  attributed  to  the  use  of  the  principle  of  linear  superposi¬ 
tion  of  airloads  in  the  flutter  analyses. 

(8)  Recently,  two-dimensional  codes  were  developed  without 
the  restrictions  of  low-frequency  (Reference  12,  19,  and  20)  and 
small-disturbance  (Reference  21).  The  viscous  effect  was  accounted 
for  by  using  a  viscous  ramp  method  (Reference  20).  Three-dimensional 
unsteady  transonic  code  has  also  been  developed  (Reference  22).  These 
new  developments  will  be  very  valuable  tools  for  exploring  the  field 
of  transonic  aeroelasticity. 

(9)  A  logical  extension  of  the  present  development  is  to 
include  a  damping  matrix  term  in  the  response  equations  of  motion  to 
simulate  the  active  control  forces  for  the  application  of  flutter 
suppression. 
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